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Abstract 

This paper revisits a homogenization problem studied by L. Tartar related to a tridi- 
mensional Stokes equation perturbed by a drift (connected to the Coriolis force). Here, a 
scalar equation and a two-dimensional Stokes equation with a L^-bounded oscillating drift 
are considered. Under higher integrability conditions the Tartar approach based on the 
oscillations test functions method applies and leads to a limit equation with an extra zero- 



< 

I order term. When the drift is only assumed to be equi-integrable in L^, the same limit 

j3 ' behaviour is obtained. However, the lack of integrability makes difficult the direct use of 

the Tartar method. A new method in the context of homogenization theory is proposed. 
It is based on a parametrix of the Laplace operator which permits to write the solution of 
the equation as a solution of a fixed point problem, and to use truncated functions even 
in the vector-valued case. On the other hand, two counter-examples which induce differ- 
Q^ ] ent homogenized zero-order terms actually show the optimality of the equi-integrability 

assumption. 

\Q . Keywords: Homogenization - Second-order elliptic equations - Drift 
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1 Introduction 

At the end of the Seventies L. Tartar developed his method based on oscillating test functions to 
deal with the homogenization of PDE's. In the particular framework of hydrodynamics [13, 14] 
he studied the Stokes equation in a bounded domain Q of M'^, perturbed by an oscillating drift 
term, i.e. 

Aus + curl (t>e) X Ue + Vps = f in 

div(M^) = infi (1.1) 
= on Q, 

where the oscillations are produced by the sequence of vector-valued functions which weakly 
converges to some v in L^{Q)^. L. Tartar proved that the limit equation of (1.1) is the Brinkman 
[4] type equation 

Am -I- curl (v) x u + Vp + Mu = f in Q. 

div(M)=0 inn (1.2) 
u = on 
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where M is a positive definite symmetric matrix-valued function. More precisely, M is defined 
by the convergences 



{Dw^fve — ^ MX weakly in L2 (0)3, for any A G M^, (1.3) 

where e W^'^{D.)^ solves the Stokes equation (1.1) in which is replaced by A. Then, the 
convergence (1.3) combined with the compactness of in L^{Q)^, yields the zero-order term 
Mu in (1.2). In [15] L. Tartar revisited this problem using the H-measures tool. On the other 
hand, the appearance of such a strange zero-order term in homogenization was also obtained 
from finely perforated domains by D. Cioranescu, F. Murat [5] for the Laplace equation, and 
by G. Allaire [2] for the Stokes equation, with zero Dirichlet boundary condition on the holes. 
Since curl {v^) x is orthogonal to u^, the energy associated with (1.1) is reduced to 

\Due\'^dx, (1.4) 

and thus does not depend on the drift Vs- Starting from this remark our aim is to study 
two drift homogenization problems associated with the same energy (1.4), and to specify the 
optimal integrability satisfied by the drift so that the Tartar approach holds. The first problem 
is scalar and the second problem is a two-dimensional equivalent of the Stokes problem (1.1). 
However, we have not succeeded in obtaining an optimal result for the three-dimensional Stokes 
equation (1.1) since the best integrability assumption for is not clear. 

In Section 2, we consider the following scalar equation in a bounded open set Q of R-^, 

— Aus + bs ■ V^e -|- div (&£ Us) = f in Q 

(1.5) 

=0 on dfl, 

where b^ G L°°(fi)^ is bounded in L'^{Q)^ . We obtain three different homogenization results: 

In Section 2.1, assuming that the divergence of the drift b^ is bounded in W~^''^{Q), with 
q > N, we prove (see Theorem 2.1) that the sequence weakly converges in Hq{Q) to the 
solution u of the equation 

Au + b ■ + div (bu) + Ltu — f in Q 

(1.6) 

u — on dQ, 

where is a nonnegative function. The proof follows the Tartar method using the oscillating 
test function 

w,:^A-'{diy{b,))eH^{n). (1.7) 

Then, in Section 2.2, assuming only the equi-integr ability of the sequence in L^(n)^ 
(this is actually a weaker assumption than the equi- integrability of the whole sequence b^), we 
obtain (see Theorem 3.1) the limit problem (1.6) with 

\Vws-Vw\'^ — ^ /i weakly in L^(ri) and fiu^ e L^{9.). (1.8) 

It seems intricate to apply directly the Tartar method with the test function Ws, since we cannot 
control the terms b^ ■ Vu^Ws and b^ ■ Vw^u^. To this end, one should consider truncations of 
both and Vw^. To overcome this difficulty we propose a new method, up to our knowledge, 
in the context of homogenization theory, based on a parametrix of the Laplace operator. It 
follows that Ue reads as a solution of a fixed point problem, which allows us to estimate the 
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sequence Vw^ ■ Vu^ only using a truncation of Vw^- The equi-integrability of Vw^ then gives 
the thesis. Also assuming that h e L^(Q)-^, with q > N, (which ensures the uniqueness in 
(1.6)) we prove the following corrector result 

Ue - {l + We-w)u — )■ strongly in WI^I{VL), for any q G [1, N'). (1.9) 

Finally, in Section 2.3, wc show the optimality of the equi-integrability condition thanks to 
a counter-example in the periodic framework (see Theorem 2.6). Making a change of functions 
with he = Vwg, equation (1.5) is shown to be equivalent to the following equation 

- AVe + lleVe ^ fe, with /I^ '-^ l^^W^f , (l-lO) 

the solution of which has the same limit as u^- G. Dal Maso, A. Garroni [6] proved that the 
class of equations of type (1.10) is stable under homogenization. Here, we do not use this 
general result, but we explicit an oscillating sequence We so that the limit equation of (1.5), or 
equivalently (1.10), is 

-Am + 7m = /, (1.11) 

with an explicit constant 7 which turns out to be < fi. Therefore, the loss of equi-integrability 
for violates the result of Section 2.2. Note that the vectorial character of the drift term in 
equation (1.5) makes difficult the derivation of a closure result similar to the one of [6] which 
is strongly based on a comparison principle. 

In Section 3, we consider the following two-dimensional equivalent of the perturbed Stokes 
problem (1.1), 

— Aue + curl (ve) Ju^ + Vpe = f in 

dw{ue) = inn (1.12) 

Ue = on Q, 

where J is the rotation matrix of angle 90°, and Ve G L°°{fl)'^ is bounded in L^(0)^. We follow 
the same scheme as in the scalar case: 

In Section 3.1, assuming that the sequence Ve is bounded in L''{Q)'^ with r > 2, we show 
(see Theorem 3.1) that the sequence Ue weakly converges in Hq{Q) to the solution u of the 
Brinkman equation 

— Am + curl (v) Ju + Vp + Mu = f in Q 

div{u) = inn (1.13) 

u — on fl, 

where M is a symmetric positive definite matrix- valued function defined by the convergence 
(1.3) inL^{nf. 

In Section 3.2, assuming only the equi-integrability of the sequence Ve in L^(f2)^, we prove 
(see Theorem 3.3) owing to the Tartar method that the sequence Ue weakly converges in i?o(0) 
to the solution u of the Brinkman equation (1.13) with similarly to (1.8), 

(Dw^fve MX weakly in L^Qf and Mu ■ u E L\n). (1.14) 

The proof is based on a double parametrix method carrying on both the velocity Ue and the 
pressure Pe- However, the proof of the last estimate of (1.14) is more delicate than the one 
of (1.8), since we cannot use a comparison principle as in the scalar case. We need to introduce 
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a test function similar to but associated with a truncation of v^. Moreover, if ft has a regular 
boundary, v e with r > 2, and M e L"*(J1)^^^ with m > 1, we get the corrector result 

u,-u-WeU — > strongly in W^'\Qf, where W^X := w^, for A G (1.15) 

Finally, in Section 3.2, we construct an oscillating sequence which is not equi-integrable 
in L^(f2)^, which leads to the limit problem (1.13) involving a matrix F which is not symmetric 
and satisfies the strict inequality 

FA • A < MA • A, for any A 7^ 0, 

which is inconsistent with the Tartar approach. This shows the optimality of the equi-integrability 
condition as in the scalar case. It would be very interesting to find the closure of the family of 
problems (1.12) under the sole condition of L^-boundedness of the sequences v^. This problem 
is far from being evident due to the absence of comparison principle for such a vector-valued 
equation. 

Notations 

2N 



The space dimension is > 2, and 2* :— 



The conjugate exponent of p > 1 is denoted by p' 
For u:R^ — > R^, Du 



N-2 

P 



p — 1 

dui 



^^J / l<i,j<N 

• For E : ^ K^x^, Div (E) ^ ^ 

V.?=i ^ / i<i<Ar 

• H^(Y), with Y := (0, 1)-^, denotes the space of the F-periodic functions on which 
belong to Hl^{M.^). 

2 A scalar equation with a drift term 

Along this section Q is a bounded regular open set of R^, with N >2, and / is a distribution 
in H-^{n). 

2.1 The classical case 

Let q G (A^, 00). Consider a sequence he in L°°(n)^ such that 

he—^h weakly in L2(0)^ and div (6,) is bounded in ^^-^'^(0). (2.1) 
Let We e Wl''^{Q) be the solution of the equation (see, e.g., [8] Theorem 2.1) 

Aw, = div(6j inD'(Q). (2.2) 
Up to a subsequence We weakly converges in Wl''\Vt) to the function w solution of 

Aw -div (6) in2)'(n). (2.3) 

We have the following result: 
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Theorem 2.1. The solution G HliVL) of the equation 

- Au^ + b, ■ Vu, + div {b^ u^) = f in D'(f^), (2.4) 
weakly converges in Hq{D,), up to a subsequence, to a solution u e Hl{Q) of the equation 

- + VM + div(6M) + = / inV'{n), (2.5) 

where /i is the function defined by the convergence 

jVu^e — Vu^l^ — ^ II weakly in hi [Q). (2.6) 

Remark 2.2. The uniqueness for equation (2.4) is not evident under tiie sole assumption 
b e L^(Q)^. Assuming a stronger integrability of b we will obtain in Theorem 2.4 the uniqueness 
for the limit equation. 

Proof. The proof is based on the Tartar method of the oscillating test functions (see Appendix 
of [12], and [16]). The function We of (2.2) will play the role of the oscillating test function. 
The variational formulation of (2.4) is 

/ Vue-V^dx+ I be-Vue(pdx- / be ■ V (fi Ue dx = {f , ip) H-^n),H^{n) , '^(peHoifl). (2.7) 
Jn Jn Jn 

Then, by the Lax-Milgram theorem there exists a unique solution Ue of (2.7) in HQ{n). In 
particular, for v e W^'°°{D,), putting (p — vu^ as test function in (2.7) we obtain the identity 

/ \Vue\'^vdx+ / Vue - Vvuedx - / b^ ■ Vv dx ^ {f , V u^) H-i{n),H^{n) , (2.8) 
Jn Jn Jn 

which will be used several times. So, choosing v = 1 in (2.8) the term with b^ cancel so that 
we easily deduce that is bounded in Hq^Q) and weakly converges, up to a subsequence, to 
a function u in Hq{Q). Therefore, it follows from (2.7) the limit variational formulation 

/ Vu-'V(pdx+ / b-'Vu(pdx+ / (pdiy- / b-Vipudx^{f,ip)H--^(^n),Hl{n), (2-9) 
Jn Jn Jn Jn 

which holds for any ip e wI''^{Vl) (due to the embedding of W^'^iVL) into C{Vl) for q > N), 
where the measure u is defined by the convergence 

be ■ Vue — ^ b-Wu + u weakly-* in M{fl). (2.10) 

The hmit equation associated with (2.9) is 

- Au + b-Wu + u + div (bu) = f in D'(Q). (2.11) 

Now, let us determine the measure z/ of (2.10). Let p G C^{Q). Putting (pw^ as test 
function in (2.7) and ipu^ in (2.2), and taking the difference of the two equalities we get 



/ ■ V</7 Wedx — / Vu^g • V</7 dx 

Jn Jn 



■ (/, (p We) H-^n),H^{n) - be- Vlie (pWedx+ / be' WWe (p Ue dx + / 6^ • V</7 Ue We dx 

Jn Jn Jn 

be • VUe (pdx — / be • V</? Ue dx. 

Jn Jn 

(2.12) 
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(2.15) 



Passing to the limit in (2.12) by using the strong convergence of in LP{fl), for p < 2*, and 
the uniform convergence of in C{Q,) {q > N), we obtain 

/ Vm • V(p w dx — Vw • V(p u dx 
Jq Jn 

— {f,^w)H-i^^^^H^(^^^ — / b-Vucpwdx— / (pwdiy+ / a(pudx+ / h-Vipuwdx 

' ° Jci Jq Jq Jq 

— / b -Vuipdx — I cpdu— / b - Vipudx, 
Jq Jq Jq 

(2.13) 

where the measure i/ is defined by (2.10) and the function a is defined, up to a subsequence, 
by the convergence 

• — ^ a weakly in L^(fi). (2-14) 
On the other hand, putting ipw E Wq''^{Q) in (2.9) and cpu E Hq{Q) in (2.3) we have 

/ Wu-V{ipw)dx = {f,(pw)H-^(n),m{n) - / b-Vuipwdx - I (pwdu 
Jq Jq Jq 

+ / b-'Vw(pudx+ / b - Vipuw dx, 
Jq Jq 

/ Vw ■V(ipu)dx ^ / b-Vu(pdx+ / b-Vipudx. (2.16) 
Jq Jq Jq 

Equating the difference between (2.15) and (2.16) to the right-hand side of (2.13), it follows 

that 

/ aipudx— / b - Vwifudx - / ipdv^Q, for any (/? e C^(Q), (2-17) 
Jq Jq Jq 

which implies that 

v^au-b-Vwu in D'(Q). (2.18) 

It thus remains to determine the limit equation (2.5). To this end, we pass to the limit by using 
if We as test function in (2.2) and the definition (2.6) of /x, and we put ipw va. (2.3), which yields 

/ {^jji + \Vw\^^ (p dx + I Vw-Vi-pwdx— i aipdx+ / b-Vcpwdx, (2-19) 
Jq Jq Jq Jq 

/ \'Vw\^ipdx+ / 'Vw-'V(pwdx= / b-'S/w(pdx+ / b-V^pwdx. (2.20) 
Jq Jq Jq Jq 

Equating (2.19) and (2.20), we deduce that 

^ = a-b-Vw inD'(O), (2.21) 

which combined with (2.18) implies that 

v^HU inD'(O). (2.22) 

Finally, the limit equation (2.11) and the relation (2.22) give the desired homogenized equa- 
tion (2.5). □ 

Remark 2.3. It can be shown that 

^x{x)= I ^iix,dOC-C, (2.23) 

where fj, denotes the matrix-valued if-measure (or micro-local defect measure) of the sequence 
be (see [15] and [7]), and S^-^ the unit sphere of . 

Assumption (2.1) is actually not sharp. In the next section we replace it by the boundedness 
of be and the equi-integrability of Vw^ in L^(Q)^. 
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2.2 The case under an equi-integrability assumption 

In this section f2 is a bounded open set of M^. Consider a sequence in L°°{fl)^ the Hodge 
decomposition of which is 

b, = Vw, + with We e H^{n), e L'^{nf and div (^e) = 0, (2.24) 

such that 

be^ b weakly in L'^{Q)^. (2.25) 

Note that for a fixed e > 0, w, E W^'P{n) and e Lp(1])^ for any p e [2, oo). But the essential 
point is the asymptotic behaviour of the sequences b^, Vw^, ^g. Our main assumption is the 
equi-integrability of the sequence Vw^ in L^(f2)^. By virtue of the Vitali-Saks theorem this is 
equivalent to the following convergence, up to an extraction of a subsequence, 

\Vw^ -Vw\^ — ^ H weakly in L^(Q), (2.26) 

(Compare to (2.6) with q > N). 
We have the following result: 

Theorem 2.4. 

i) Under the equi-integrability assumption (2.26) the solution of (2.4) weakly converges in 
Hq^Q) to a solution u of the equation 

- Au + b-Vu + div {bu) + fxu^ f inT)'{Q), (2.27) 

with 

/ l^u^dx<{f,u)H-i(n),H^{n)- / \^u\^ dx. (2.28) 
Jn Jn 

ii) Also assume that b e L^(Q)^, where q > 2 if N — 2 and q — N if N > 2. Then, we have 

/ \Vu\'^dx+ / f^u^dx^ {f,u)H-i^a),m{n)- (2-29) 
Jn Jn 

and there exists a unique solution u G Hq{Q) of equation (2.27), with jiu^ E L^iVl). 
Moreover, for any p e [1,2) if N = 2 and p = N' if N > 2, we have the corrector result 

Vue -Vu- {Vwe - Vw) u — > Q strongly in Lf„^(Q)^, (2.30) 

and for any r e [l,p), 

Ue-{l + We-w)u — )■ strongly in W^^^{n). (2.31) 

Remark 2.5. No equi-integrability is required for the divergence free sequence ^g. Actually, we 
can prove that the equi-integrability of the sequence b^ in L^(fi)^ implies the equi-integrability 
of its two components Vw^, in ^^^{fl)^. Therefore, condition (2.26) is really weaker than 
the equi-integrability of b^. 

Moreover, the equi-integrability of Vw^ in L^(f2)^ is essential for deriving the limit equation 
with the zero-order term /i u. When this condition is not satisfied we can obtain a similar limit 
equation but with a different zero-order term (see Section 2.3). 
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Proof of Theorem 2.4. The limit u of Ue in HliVt) solves the equation (2.11) where v is 
defined by 

be ■ Viie -h-Vu — ^ V weakly-* in M(Q), (2.32) 
By the Murat, Tartar div-curl lemma [10] we have 

he ■ Vue = i^e + Vw) ■ Vue + (Vw^ - Vw) ■ V Ue — ^ h-Vu^v in T)'{Q). 

This combined with the equi-integrability of Vw^ implies that v is also given by the convergence 

{Vwe - Vw) ■ Vue — ^ y weakly in L^{Vt). (2.33) 



The proof of Theorem 2.4. is based on a parametrix method which allows us to express Ue 
as a solution of a fixed point problem. As a consequence, we obtain a strong estimate of Vue 
in L\^^{VL) for some p > 1 close to 1. However, this estimate cannot provide directly the desired 
limit V of (2.33) since p < 2. To overcome this difficulty we consider a truncation t]^ of Vw^ 
which is bounded by A; > 0. Then, we can pass to the limit as e tends to zero in the product 
T]^ ■ Vue for a fixed k. Hence, thanks to the equi-integrability of Vwe we deduce the fimit v as 
k tends to infinity. 

The proof is divided into four steps. In the first step we present the parametrix method 
which leads to a L^-strong estimate of Vm^. In the second step we determine the limit of the 
sequence r}^ ■ Vm^ for a fixed k > 0. In the third step we determine the limit and the limit 
equation (2.27) together with (2.28). The fourth step is devoted to the proof of equality (2.29) 
and the corrector results (2.30) and (2.31). 

First step: The parametrix method. 

First, let us define a parametrix for the Laplace operator in Q. To this end consider two 
sequences of functions ^Pmi^n in C'^(Q), such that 

< v^n, "^n < 1 and = 1 in supp (■?/'„) , for any n > 1, 

[n>l: supp (-^n) n X 7^ 0} is finite, for any compact subset K <zVL, , . 

[2.64:) 

ipn — I in Q. 

n>l 

Let E be the fundamental solution of the Laplace operator in R^. Then, the operator P defined 
in D'(n) by 

^C:=$^^n£^*(<^nC), forCeD'(O), (2.35) 

n>l 

is a parametrix of the Laplace operator (see [1] Chapter I, for further details) which satisfies 
P{AC)^C-KC and A(PC) = C - ^^'C, for C e 2)'(n), (2.36) 

where K, K' are two C°°-kernel operators properly supported in Q. Thanks to the Calderon- 
Zygmund regularity for the Laplace operator (see, e.g., [8] Theorem 2.1, and the references 
therein) we also have for any p > 1, and s e [0, 2] such that s -|- ^ is not an integer, 

P maps continuously D'(l]) to D'(l]), and W^^^'^i^) to Wl^''''{n). (2.37) 
Then, applying (2.36) to the solution Ue of (2.4) we have 

Ue = P{AUe) + KUe = P [div(w {VWe - Vw) )] + P [div(Vw£ {Ue - u))] 

(2.38) 

+ P [div {uVw) ]+P{ie- yUe + he ■ WUe - /) + KUe, 
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Fix p > 1 close enough to 1 and s e {N/p', 1). Since — u strongly converges to in L'^{Q) 
for any q G (2, 2*), the sequence div(Vw£ {u^ — u)) strongly converges to in W~^'P{Q), hence 
by (2.37) we have 

P [div(Vwe (ue-u))] — >0 Strongly in W^^'^{n). 

Moreover, the sequence ■ Vm^ + 6^ ■ Vm^ is bounded in L^{Q), thus in W^''^'^{Q) since s > N/p'. 
Therefore, again by (2.37) the sequence VP(^£ ■ Vm^ + ■ Vu^ — /) is bounded in W^~^'P{fl)'^ , 
and up to a subsequence strongly converges in L^^^{fl)^ . Hence, since 

■ + he ■ Vu^ — ^ ^ • Vli + + 6 • Vii in D'(0), 

we deduce from (2.38) the strong estimate 

VUe - VP [div(M {VWe - V w) )] 

= VP [div {uVw) + ^-Vu + v + h-Vu- f]+ VKu + OiP^(f2)^(l) (2.39) 
= VP[v + h- Vu + d\Y{hu) - /] + V-fr'u + OiP^(n)iv(l) (since ^-Vu^ d\Y{u^)), 

where Oiv (j-j)iv(l) denotes a sequence which strongly converges to in L\^^{yL)^ . On the other 
hand, by (2.36) and (2.37) we have 

VP [div(M {Vwe - Vw) )] = VP [A(m {we-w))]- VP [div(VM (w^ - w))] 

= VP [^{u {w, -w))]+ OiP^(n)^(l) 

= V(li (We - W) ) + OiP^(Q)Jv(l) 
^u{VWe - Vw) + OiP^^(Q)iv(l). 

Therefore, this combined with (2.39) yields 

Vue - u {Vwe - Vw) = VP [i/ + 6 • Vii + div {hu) - /] + VKu + olp^^^q)n{1). (2.40) 

Second step: Estimate of the sequence 77^ • Vu^. 

Set Tj^ := Vwe l{|vu;e|<fc}, for a positive integer k. Let us determine the limit of 77^ • Vu^ in 
-^k)c(^)^- Using a diagonal extraction, there exists to a subsequence of e, still denoted by e, such 
that 7]'^ weakly converges to some 7]^ in L°°(r2)^ for any k. By the strong convergence (2.40) 
combined with the weak convergence of u (Vw^ — Vw) to in L^iVL)^ (for p close to 1) we have 

77^ • Vwe - (?7^ - jy'^) • (Vwe - Vty) w 
— ^ T]^ ■VP[v + h-Vu + div (6ii) - /] + 77*^ 

Hence, we get that 

(7^ ^ u + r]'' ■ VP [i/ + 6 • Vii + div (6ii) 

where 

/i'' := lira [(r?^ - r?'') ■ (Vw^ - Vw)] 
Third step: Determination of v and the limit equation (2.27). 



VKu weakly in L\^^{Q). 

- /] + 77^^ • VKu in Q, (2.41) 
weakly in -L^(il), 

(2.42) 

weakly in L'^iVL). 
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Starting from the limit equation (2.11) we have by (2.36) 

P[u + b-Vu + div{bu) - f]+Ku in n, 

hence 

rj'' ■Vu^r]'' -VPlu + b-Vu + diyibu) - f] ^-q^-VKu in Q. 
Equating this with (2.41) we obtain 

(T^ ^ f/u + f]^ -Vu inO. (2.43) 

Now, let us pass to the limit as A; — > cxd. By virtue of the equi-integrability of Vw^ in L?{yt)^ 
and by definition (2.42) the sequence strongly converges in V-{VL) to the function of (2.26), 
Tj^ strongly converges to Vw in L^(r2)^, and strongly converges to v + Vw ■ Vu in L^(f2). 
Then, up to a subsequence converges to a.e. in Q, and by the Fatou lemma combined 
with equality (2.43) we get 

/ < liminf / l/x*^ m| (ix < liminf / — ■Vu\dx — I \v\dx. (2.44) 

Jq Jn '^-^"o Jn Jn 

We deduce from (2.44) and (2.43) that i^ue L^{fl) and 

u^fjLU in Q, (2-45) 

which yields the limit equation (2.27). 

It remains to prove the inequality of (2.28). Let v e L°°(Q) and t e R. By (2.26), (2.33) 
and (2.45) we have 



— Vu — ( — Vw) tv)\^ dx 



I 

Jn 

= / \Vue -Vul"^ dx + f / \Vw^ -Vw\'^v^ dx -2t / Vue ■ {Vwe - Vw) v dx + o{l) 
Jn Jn Jn 

= - / \Vu\^dx + t'^ / nv'^dx-2t I iJ,uvdx + o{l), 

Jn Jn Jn 



(2.46) 



hence 



/ iJ,v^dx — 2t / iJiuv dx -\- {f ,u) H-^{n),H^{n) ~ / | Vxi|^ > 0, e R. 
Jn Jn ' ° Jn 

This implies that 

Let Tfc, A; > 0, be a function in C^(M) such that 

o<ii<i and ;[|*|i'+2^ 

Putting V — Tk{u) as test function in (2.47) and using that Tk{uf' < uTk{u), we get 
liuTk{u)dx^ < (^{f,u)H-^Q),H^{n) - J \^u\'^dx^ J iiuTk{u)dx, 
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hence 

/ fiuTkiu)dx<{f,u)H-i(n),Hl^{a)- / \^u\^ dx. (2.49) 
Jn Jn 

Since uTk{u) is a nondecreasing nonnegative sequence which converges to a.e. in the 

Beppo-Levi theorem apphed to (2.49) thus gives inequahty (2.28). 

Fourth step: Proof of equahty (2.29) and of the corrector results (2.30), (2.31). 

Assume that h G L«(fi)^, where g > 2 if = 2 and g = iV if iV > 2. Let be a sequence in 
Co(M) which strongly converges to u in i^o(fi) and a.e. in fi, and such that \Vipn\ is dominated 
by a fixed function in L^(Q). Putting the truncation function Tf;{(p„) (2.48) in the limit equation 
(2.27) we have 

/ Vu ■ VTk{(pn) dx + / b-VuTk{ipn)dx - / b ■ VTk{(pn} u dx + / jjLuTk{ipn) dx 
Jn Jn Jn Jn 

Since b • Vm, ij,u G L^(fi) and bu G L^(fi)^ (as a consequence of 6 G L'^{n)^), we can pass 
to the limit as n — > oo in the previous equality owing to the Lebesgue dominated convergence 
theorem, which yields 



\/u-\/Tk{u)dx+ / b-VuTk{u)dx— / b-VTk{u)udx+ / iiuTk{u)dx 

Jn Jn Jn (2.50) 

= (/ , ^fc('"'))if-i(f^),ifi(n)- 

Then, using that |T'fc(M)| < \u\, < T'^(m) < 1, Tk{u) strongly converges to u in Hq^Q), and 
that felt G L^(Q)'^, fjLU^ e -L^(Jl), and passing to the limit as /c — >■ oo owing to the Lebesgue 
dominated convergence theorem we get 



/ \\/u\'^ dx + / b-Vuudx— I b -Vnudx + / fiu^ dx = {f,u) H-^{n),H^ 
Jn Jn Jn Jn ' ° 



(n) , 



which is (2.29). Moreover, the proof of equality (2.29) with / = shows that there exists a 
unique solution u G Hq^Q) of equation (2.27), with fiu^ E L^{Q). 

It remains to prove the corrector results. By the estimate (2.46) with v = T^iu) and t = 1, 
combined with equality (2.29) we have 



^lim lira (^j \Vu^ — Vu — (Vu^^ — Vw) Tk{u)\^ dx^ 

— lim ( / jiiu — Tdu)Y dx \ — 
''-^'^ \Jn J 



(2.51) 



On the other hand, let p G [1, 2) if = 2 and p — N' \l N > 2, and consider an open set 
a; (S Q. By the Holder inequality we have 

/ I VWe — Vm — {WWe — Vw) dx 

J U) 

<2P~M /" \Vu,-Vu- {Vwe-Vw)Tkiu)f' + [ \Vw,-Vwf \u-Tk{u)fdx] 
\Jn Jul / 

p p 

<c^y" \Vue-Vu- {Vwe-Vw)Tk{u)\'^Y +^(^J \u - Tk{u)\^ dx^ 



(2.52) 



£ , ^ 1 _E 

2 / r 2p \ 2 



<c[ I \Vue-Vu- {Vwe-Vw)Tk{u)\ +c( / lul'^dx 

/ \J{\u\>k}nu> 
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Since u e L2-p(a;) by the Sobolev embedding, passing successively to the hmits e ^ and 
/c — > oo in (2.52) owing to convergence (2.51) we obtain the strong convergence (2.30). 

Let r e [l,p)- Since w^ — w strongly converges to in L^(uj), by the Holder inequahty the 
sequence {w^ — w)'Vu strongly converges to in ^{ijj)^. Finally, this combined with (2.30) 
implies the corrector result (2.31). 



2.3 A counter-example 

In this section 17 is a regular bounded open set of R^, and Y :— (— |, |)^. For fixed R G (0, |) 
and A* > 0, let e (0, R) be defined by the equality 

= II. (2.53) 



|lnr£| 

Let Ws be the F-periodic function and be the sF-periodic function defined by 
' In r — In 

- n r := \y\ t yr^.n) 

yeY, We(x) := (^^) , xE R^. (2.54) 



lni?-lnr. ^ :^ \y\ e {r,, R) 



e 

si r < r, 



1 si r > R, 

Note that by (2.53) we have 



We then consider the drift be defined by 

b,{x) := Vw,{x) = ^ VWe , for x e (2.56) 

Taking into account (2.53) it is easy to check that 

We — ^ 1 weakly in H^{il) and weakly-* in L°°{il). (2.57) 

Let / be a non-zero function in L^(Q). We study the asymptotic behavior of the equation 
(2.4) with the drift be of (2.56), i.e. 

-Aue + Wwe-Vue + div{VweUe)^f m V (fl) . (2.58) 

We have the following result: 

Theorem 2.6. The solution Ue of (2.58) weakly converges in Hl{Q) to the solution u of the 
equation 

3 ie^ - I) 

-Au + -fu^f in'D'in), where 7:=—^- u < u. (2.59) 

V ' 4 (e2 -I- 1) ' ^ ^ 

Remark 2.7. Using the periodicity we can check that the sequence l^^p = |Vtfep converges 
in the weak-* sense of measures on fl - but not weakly in L^{fl) - to the constant /j, defined by 
(2.53). Theorem 2.6 can thus be regarded as a counter-example to the statement of Theorem 2.4 
without the equi-integrability assumption on the drift be in L^(Q)^. Indeed, the conclusion of 
Theorem 2.4 would give a hmit equation (2.59), with 7 = /x. 
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Proof of Theorem 2.6. The proof is divided into two steps. In the first step we construct an 
oscillating test function which solves equation (2.64) below. In the second step we determine 
the limit equation (2.59). 

First step: Construction of an oscillating test function. 

Denote by Qr the disk of radius r centered at the origin. Consider the unique solution in 
H^{Qr) of the equation 

-^Az, + i|viy,|'z, ing^ 



dZ, 



dn 



(2.60) 

= on dQji. 



The function Z^ is radial and can be computed explicitly. Using the Laplace operator in polar 
coordinates and iVVFep = a^r"^ 1q^\q^^, we get 



^ „2 



r^ + Ce if r e (0, r^] 



Ze{r) = { 2 where , , y (2.61) 

The constants Og, fog, Cg are determined owing to the boundary condition on BQr and to the 
transmission conditions on dQr^, i.e. 

Z'^{R)^Q and Z,(r+) = Z,(r,-), Z^(r+) = Z^(r,-). (2.62) 

We extend Z^ by the constant value Z^{R) in F \ Qr, and by y-periodicity in the whole 
space R^. The y-periodic extension is still denoted by Z^. An explicit computation combined 
with (2.53) yields 

Zg^Z:=i|^^±^i strongly in i/;(y). (2.63) 
As a consequence of (2.60), (2.61) the rescaled function Ze{x) := Z^[^) is solution of the equation 

-Az, + |V^/;,rZg = xL(7) m2)W, (2.64) 

where Xg^ is the F-periodic function agreeing with in the period cell Y . Moreover, the 
following convergences hold 

z^^ Z weakly in H^{9) and Xq^ (^) ^ 1 weakly-* in L°°(fi), (2.65) 

where the constant Z is defined by (2.63). 

Second step: Determination of the limit equation (2.59). 

Define the function := e^~^^ u^. Then, equation (2.58) is equivalent to 

-l^Ve + \Vwe\^Ve^e^-'"' f inD'(Q). (2.66) 

G. Dal Maso, A. Garroni [6] proved that this class of equations is stable under homogenization. 
In the present case, the use of the oscillating test function z^ will allow us to obtain the limit 
equation (2.59). 
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On the one hand, choosing i> = in (2.8) we get 

/ \'S7We\'^ ul dx — / -VUsUedx — / \'S/Ue\'^ dx — / fWeUedx<C, (2.67) 

Jn Jn Jq Jn 

since is bounded in Hq{D,) and < < 1- Then, by the Cauchy-Schwarz inequahty we 
have 



Jn \Jq, ) \Jn 



VwJ'^u'tdx <c + c( / |Vii£pdx) ( / iVwJ^utdx 



(2.68) 



<c+c' i I \Vws\^u^^dx ' 

n 



hence MgVwg is bounded in L^(0)^. This combined with convergence (2.57) imphes that 
weakly converges to u in Hq^Q). 

On the other hand, for cp e C^{il), putting the functions cp in (2.66) and ^pv^ in (2.64), 
taking the difference of the two equahties, and passing to the hmit owing to convergences (2.65) 
we obtain the equahty 

Vu-VipZdx+ ipudx^ fipZdx, for any e C~(0). (2.69) 
Jn Jn Jn 

which is the variational formulation of equation (2.59), with 7 = Z~^. □ 

3 A Stokes equation with a drift term 

3.1 The classical case 

In [13, 14] L. Tartar noted that the nonlinear term of the three-dimensional Navier-Stokes 
equation for the divergence free velocity u reads as 

(u-V)u^ Div (u^u) ^ curl (u) x u + V \u\'^) . (3.1) 

This led him to study the perturbed Stokes equation 

-Au + curl (v) xu + Vp^ /, (3.2) 

where a given vector-valued function v replaced the velocity u of the Navier-Stokes equation. 
The equivalent of transformation (3.1) in two-dimension is 

Div {u^u) = curl {u) -|- V (| , 

^0 -1\ (3.3) 



where cut\(u) :— d\U2 — d2U\ and J := 



1 



More generally equality (3.3) extends for any divergence free functions v to the following one 

curl [v) Ju = Div (u ® m) + {Duf' v — V {v ■ u) . (3.4) 
Similarly to (3.2) this leads us to the two-dimensional perturbed Stokes equation 

- Aii curl {v) Ju + Vp = f. (3.5) 
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Let be a bounded domain of R^. Let Vg be a sequence in L°°(f2)^ and let / be a distribution 
in H~^{Qy. Consider the perturbed Stokes equation 

- Aue + curl (vs) Jue + Vps = f in f2 

div (ue) =0 in O (3.6) 
Ue — on dQ. 

In the three-dimensional case where curl (v^) xu^ replaces curl (v^) Ju^, L. Tartar [14] derived 

a Stokes equation with a Brinkman law under the assumption that is bounded in L^{Q)^ 
(see Introduction). Mimicking the Tartar approach in dimension two we can derive a similar 
homogenized equation using the test function w^, for A e M^, solution of the Stokes equation 

-Aw^ + Dw{(ve-v)®X)+Vq^ =0 in Q 

div(w^) =0 inQ (3.7) 

=0 on on. 

Then, we have the following result: 

Theorem 3.1. Assume that is bounded in L^{Q)^, with r > 2. Then, the solution of 
(3.6) weakly converges in HliVL) to the solution u of the Brinkman equation 

Au + curl (v) Ju + Vp + Mu = / inQ, 

div(ii) =0 inn (3.8) 
M = on dfl, 

where M is the positive definite symmetric matrix-valued function defined by 

{Dw^Yv^ — ^ MA weakly in L^(ny and in Lf (n^ 

, /orA,/xeM. (3.9) 

Dw^ ■ Dw^ — ^ MX- n weakly-* m MiVlf and m L^^iS^f , 

Moreover, the zero-order term of (3.8) is given by the convergences 

{Dus)'^{vs ~ v) — ^ Mu weakly in {^ly 

^ (3.10) 
Due : Dw^ — ^ Mu-X weakly-* in M(0) and in L^^ {nf. 

Proof. By the representation formula (3.4) we have 

curl {Ve) JUe = {DUef + Div (u^ (g) Me) - V {v^ ■ U^) ■ (3.11) 

Hence, the variational formulation of (3.6) reads as 

/ Due : Dipdx+ / {Due fve -(fidx- / {v, (g) u^) : Dipdx = (/, u^) H-^a)^,H^{a)2 , 
JO. Jn Jn 

for any (p e H^{9.y, div (</?) = 0. 



(3.12) 



By the Lax-Milgram theorem there exists a unique divergence free function Ue e Hq{Q.)'^ solu- 
tion of (3.12). Then, putting the velocity as test function in (3.12) it follows that 



/ 

Jq 



\Due\'^ dx = (/, Ue)H-^{nY,Hl{nY , (3.13) 
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which imphes that is bounded in HQ^fiy. Let a; be a regular domain of fl. Applying (3.12) 
to divergence free functions in Hq{ujY, there exists a unique in L^(a;)/]R such that equation 
(3.6) holds in 2)'(a;)^. Moreover, by (3.11) and the boundedness of in L^(Q)^ the sequence 
Vps is bounded in H^^{uY. Hence, due to the regularity of u the sequence is bounded in 
L'^ioj). Then, considering an exhaustive sequence of regular domains the union of which is ^2, we 
can construct in a pressure p^ which is bounded in Lf^^iVL). Therefore, up to a subsequence 
the following convergences hold 

Up — ^ u weakly in HliVLY 

° (3.14) 
Ps — ^p weakly in L2^Qj,(f])/M, 

Now, in view of (3.11) it is enough to determine the limit of the term (Due)^ v^- By the 
regularity results for the Stokes equation (see, e.g., [9] Theorem 2, p. 67) the sequences and 
satisfy 



weakly in H\n)^ and in W^/^iQf 

(3.15) 

^ weakly in L2(0)/K and in L[,^(0)/M. 

which imply convergence (3.9). Let ip G C^{Q). Following the Tartar method we put (fiw^ in 
equation (3.6) and Lpu^ in equation (3.7). Then, from the representation (3.11), the conver- 
gences (3.14), (3.15) and the boundedness of in L^{fl) we deduce that 

/ Due : Dw^ ifdx- / (v^ (g) u^) : Dw^ dx = o(l) 

J? ^ for any ip E C^in), 

/ Dw^ : Du^ (p dx — / (^{v^ — v) <Si X) : Dus(p dx — o{l), 
Jn Jn 



hence 



A f r)/>,.A^T„^^ _ .^^ 

in D'(Q). (3.16) 



Due : Dw^ - {Dw^ fv, ■ u^ — ^ 



{Du,f{v, -v)-X- {Dw^fv, -u, ^0 

By virtue of the strong convergence of u^ in any L*(fi)^ space for s G (1, oo). convergences (3.16) 
and (3.9) imply (3.10). This combined with (3.11) yields finally the limit problem (3.8). □ 

Remark 3.2. It can be shown that 

M{x) = [ [tr (/. (x, dO) - M (x, e • e] e ® e (3.17) 

where /it is the matrix- valued iJ- measure of the sequence (see [15, 16]). 

The case where is only bounded in L^(r2)^ is much more delicate. On the one hand, 
under additional assumptions we will extend the Tartar result when is bounded and equi- 
integrable in L'^{fiy. On the other hand, we will give an example of a sequence for which 
the homogenized Brinkman equation is not the one obtained by the Tartar procedure. 

3.2 The case under an equi-integrability condition 

In this section we make the following weaker assumption on the drift, 

Vs — ^ V weakly in L^(f2)^ and is equi-integrable in L^(l])^. (3.18) 
Then, we have the following extension of Theorem 3.1: 
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Theorem 3.3. 

i) Under the equi-integrability assumption (3.18) the solution of (2.4) weakly converges in 
Hq{Q) to the solution u of equation (3.8) with 

/ Mu-udx<{f,u)H-i{Q)2m^n)2- / \Du\'^dx, (3.19) 
Jn .In 

where M is the positive definite symmetric matrix-valued function defined by 
(Dw^fv, MX weakly in LVQ)^ 

\ ^ ^ forX,iieR^ (3.20) 

Dw^ : Dw^ — ^ MX- 11 weakly-* in M(Q)^ 

ii) Also assume that Vt has a Lipschitz boundary, v e U{fiy, with r > 2, and M e L"^{fiy^'^, 
with m > 1. Then, we have the equality 

/ \Du\'^ dx+ / Mu-udx^ (/, u)H-i{Q)2,H^(ny , (3-21) 

and there exists a unique solution u G -ffo(^)^ of equation (3.8), with Mu ■ u £ L^iVt). 
Moreover, we have the corrector result 

Ue-u-WeU — 5^0 strongly in W^'^{Vtf, (3.22) 

where is the matrix-valued function defined by 

W,X:^w^, forXeM.^. (3.23) 

Remark 3.4. Contrary to Theorem 2.4, in the part ii) of Theorem 3.3 we need to assume a 
higher integr ability for the matrix- valued M. Indeed, we cannot apply a truncation principle 
on Mu ■ u. Moreover, the regularity of fl is necessary to obtain the density of the smooth 
divergence free functions in the space of the divergence free functions of Hq{Q,)'^. 

Proof of Theorem 3.3. As in the proof of Theorem 3.3 the sequence is bounded in ifg (Q)^, 
and thus in any L^(f2)^ space. Then, in view of (3.11) and (3.6) together with the boundedness 
of and Vs the sequence Vps is bounded in L^(f2)^ + W~'^'^{Qy for any r G (1,2). Hence, 
thanks to the embedding of Ll^^{fl) into W^^'^{fl) for any r > 1 and a > 2/r', the sequence 
is bounded in L[qp(0)/M for any r G (1,2). Therefore, up to a subsequence we have the 
convergences 

Us — ^ u weakly in H^(0,)'^ 

°^ ^ (3.24) 
Ps — ^ p weakly in L'[^^{ft)/R, for any r G (1,2). 

The problem is to determine the vector- valued distribution i/ defined by 

curl (v^) Ju, - curl (v) Ju ^ u in D^nf. (3.25) 

Taking into account the representation formula (3.11) and the equi-integrability of in L^(Q)^, 
u is actually in L^(0)^, and is given by 

{Du^five -v)^ V weakly in L^(Q)^ (3.26) 

so that u is solution of the equation 

- Aii + i/ + curl {v)Ju + Vp^f in D'(Q) . (3.27) 
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Prom now on the proof follows the same scheme as the one of Theorem 2.4 using a repre- 
sentation of the velocity and the pressure owing to the parametrix P of (2.35). The proof is 
divided into fifth steps. The first step deals with a double parametrix method for both and 
Pg, which allows us to derive a strong approximation of Du;,. In the second step we compute 
the limit of the sequence {Du^)^v^, where is a truncation of for a fixed k > 0. In the 
third step we obtain the limit equation (3.8). In the fourth step we prove inequality (3.19). 
The fifth step is devoted to the proof of equahty (3.19) and the corrector result (3.22). 

First step: The double parametrix method. 

Consider the parametrix P (2.35) for the Laplace operator. Abusively we denote by A the 
vector-valued Laplace operator as well as by P the associated vector-valued parametrix each 
component of which is defined by (2.35). Taking the divergence of equation (3.6) we have 

= div (/) — div (curl (v^) Ju^) in 

hence by (2.36) 

Pe^P [div (/) - div (curl {v^) Ju^) ] + Kp, in Q. (3.28) 
Substituting p^ by the right-hand side of (3.28) in (3.6) it follows that 

= curl (v,) Ju, - VP (div (curl (v,) Ju,) ) + VP (div (/))-/ + VKp, in Q, 

hence again by (2.36) we have in 

u,^P [curl (ve) Jue - VP (div(curl {Ve) Jue))] + P[VP(div (/) ) - /] + L(m„J9,), (3.29) 

where L is a C°°-kernel operator acting on the pair [ueiPe)- Using the representation (3.11) of 
curl {ve) Jue-i and setting 

:= Div((ve -v)® Ue) - V{{v, - v) ■ Ue), (3.30) 

we get 

u,^P [{Du.fv, + ge-VP {diy{{Du,fv, + g,))] + F{u„p,), (3.31) 

where 

P(C, 9) := P [Div(t; ® () - V(^; • C) - / - VP (div(Div(^; ® C) - V(^; • C) - /))] 

+ L{c,e). 

Note that by (2.37) we have 



(3.32) 



Fius^Pe) — > F{u,p) strongly in W^^^{^), for any r G (1,2). 

Moreover, by (3.26) the sequence {Du^)'^Vg weakly converges to z/ + {Du)'^v in L^(r2)^ which is 
compactly embedded in W^^'^{^y for any r G (1,2). Hence, as in the first step of the proof of 
Theorem 2.4, from (3.31) and the two previous convergences we deduce, for any r G (1, 2), the 
strong convergence 

u,-P [g, - VP (div (ge) )] P [ly + {Dufv - VP (div(i/ + (Dufv))] + F{u,p) 

strongly in Wl^liVLf. 



(3.33) 



Second step: Determination of the limit a'^ of {DueYv 



T„,k 
e ■ 
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Fix r e (1,2) such that (3.33) holds. Set 

z,:^P[g,-VP{dw{g,))] and := P(div (5,) ) . (3.34) 
In view of (3.30) the sequence weakly converges to in W~^''^{0,y, hence by (2.37) we have 

weakly in W^^^{nf and Qe ^ weakly in L[„^(f])/M. (3.35) 
Moreover, by (2.36) we have 

Aze = ge- Vqe - K'ge and Ag^ = div {ge) - K'q^ in Q, (3.36) 

hence 

A (div (ze) ) = K'q, - div (K'g,) strongly in L[„,(l^)', say. 
This combined with the first convergence of (3.35) and (2.37) yields 

div (z,) — > strongly in W^^^{nf. (3.37) 

On the other hand, set := '^{\v^\<k}-, foi' a positive integer k. Up to a subsequence of e 

still denoted by £, weakly converges to some function in L^(r2)^ for any k. Consider for 
A e IR^, the solutions w^'^ and g^'*^ of the Stokes problem 



- Aw^''^ + Div {{v'l - v'') ®X)+ Vg,^'^' = in O 

div {w^'^) = in Q (3.38) 

which consists in a approximation of equation (3.7). By the regularity results for the Stokes 
equation (see, e.g., [9]) we have 

w^^^ weakly in <;^(fi)2 

for any s e (1, 00). (3.39) 

qKk ^0 weakly in Lf„^(n)/]R, 

Choose s :— r', and apply the Tartar method (see Appendix of [12]). Let 9? G C^{Vi). Putting 
(fiw^''' in the first equation of (3.36) and (pz^ in equation (3.38), and using the definition (3.30) 
of ge and the convergences (3.35), (3.37), (3.39) we have 

{Dz,f{v^^ - v^) • A - {Dw^'''f{ve -v)-Ue-^ weakly in 2)'(fi)l 

Hence, since Dz^ weakly converges to in L''(n)^^^, we deduce that 

{Dzefv^ M^u weakly in D'(0)^ (3.40) 

where the matrix- valued function is defined by 

[Dw^'^fv, M^X weakly in Lf„^(Q)^ for any s e [1, 2). (3.41) 

Now, we are able to determine the limit of the sequence {Duef'V^ in L^(0)^. With the 
definition (3.34) of z^ the strong convergence (3.33) implies that 

{Duefv^^ - (Dzefv^ {DP [u + {Dufv - VP (div(i/ + {DuYv))]f + {DF{u,p)fv^ 

weakly in L\^^{Q)'^. 
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This combined with (3.40) thus yields 

= M^u+{DP [v+{DuYv - VP {diw {v + {Duf v))]Y + {DF{u,p)fv'' . (3.42) 

Third step: Determination of the hmit equation (3.8). 

The function u solves the equation (3.27) which by (3.4) and similarly to (3.31), can read as 
u^P[u+ {Dufv - VP (div(i/ + {Dufv))] + F{u,p). 

This imphes that 

(Dm)V = {DP [v+{Dufv - VP {<liw{v + {Dufv))]f Vk + {DF{u,p))^v\ 

Therefore, equating the previous equation with (3.42) yields 

a'' = (Dufv'' + M^u in Q. (3.43) 

It remains to pass to the limit as k tends to infinity. Due to the equi-integrability of in 
L^(Q)^ and by convergence (3.26) the sequence a'^ strongly converges io v + {Du)'^v in L^(f2). 
On the other hand, putting the function w^''^ — both in equations (3.7) and (3.38) we get 
the equality 

/ \Dw^^^ - Dw^l'^ dx = [ {Dw^/ -Dw^f{vl-Ve-v'' + v)-\dx, 
Jn Jn 

which, again by the equi-integrability of v^, yields 

lim sup ( / \Dw^''' - Dw^l'^ dx] = 0. (3.44) 

Estimate (3.44) implies that the sequence defined by (3.41) strongly converges in L^[n)^^^ 
to the matrix- valued function M defined by (3.20). In particular, up to a subsequence M'^ 
converges to M a.e. in Q. Then, by the Fatou lemma combined with (3.43) and the strong 
convergences of a'' in L^(Q)^ and Vk in L^(Q)^, we get that the function Mu belongs to L^{QY. 
Finally, passing to the limit in (3.43) we obtain the equality 

u — Mu in Q, 

which gives the limit equation (3.8). 
Fourth step: Proof of inequality (3.19). 

Similarly to (3.23) let W^, k > 0, he the matrix-valued function defined by := w^''', 

where w^'^ solves (3.38). We simply denote w*'^ when A = Cj := (2 — i, i — 1), for i = 1, 2. Let 
if e Cl{nf, and let t e M. Using (3.13) we have 



/ \Du, -Du-t D{W^ip)\^ dx = (/, u)H-i(n)^H^{n)^ - / \Du\^ dx 
Jn Jn 

-2t [ Du,:D{W^cp)dx + t'^ I \D{W^^p)f dx + o{l). 
Jn Jn 



(3.45) 



Moreover, similarly to the second convergences of (3.9) and (3.10), we have for i,j = 1,2, 
Du^ : Dwp — ^ M'^u ■ Ci 

weakly in Lf„,(f]), for any s G [1, 2), (3.46) 

Dw'/ ■ Dwl'" — ^ M'^Ci ■ e-j, 
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where (compare to the definition (3.41) of M*^) the matrix- valued M'^ is defined by 

{Dwi'''fv^ ^ M'^ei weakly in (0)2. (3.47) 
Then, from convergences (3.39) and (3.46) we deduce that 

2 

/ Due : D{W^ip) dx [ Du, : Dwf ifi dx + o(l) — > I M^u ■ dx, 

2 

/ \D{W^Lp)?dx = V / Dw'/ : Dwi'^ipiipjdx + o{l) — y [ M'^ip-ipdx, 

This combined with (3.45) implies that 

/ \Du^-Du-tD{yV^ip)^dx = l^f,u)^-^(^^-)2^Hii^^-)2- \ \Du\^dx 
Jn Jci 

-2t M''u-^dx + t^ / M^<^- ^dx + o{l). 
Jn Jn 

Therefore, we have for any i e R, 

(p-(pdx-2t / M^K • (^dx + (/, 1^)^-1(^)2 ^^i(Q)2 - / iDul"^ dx >0, 
Jn Jn ' Jn 



(3.48) 



hence 

2 



Let 5 > 0, and let u be an open set such that u ^ fl. Since by (3.41) and (3.47) M'^ and 
M'^ belong to L^{ujy^'^ for s e [1,2), putting in (3.49) strong approximations cp of j^^^ in 
L2^'(Q)2, we get 

< [^{f,u)H-^nr,H^^inr- JjDufdxj j^-^^—^^dx, 



which by the arbitrariness of cu yields the inequality 

2 



Recall that, by virtue of the equi-integrability of in L'^(il)'^, the sequences M'' and 
strongly converge to M in L^{Q)'^^'^, thus converge, up to a subsequence of /c, a.e. in Q and 
in a dominated way. Therefore, passing to the limit as /c — > oo owing to the Fatou lemma for 
the left-hand side of (3.50) and owing to the Lebesgue dominated convergence theorem for the 
right-hand side of (3.50), it follows that 
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which imphes the inequality 

/ i^^r 1^1 < (/,-») g-i(n)2,ifi(n)2 - / \Du\^dx. (3.51) 

Jq \U\ Jq 

Finally, applying the Fatou lemma in (3.51) as (5 ^ we obtain the desired inequality (3.19). 
Fifth step: Proof of equality (3.21) and of the corrector result (3.22). 

Assume that Q has a Lipschitz boundary, v G L^{Q)^ , with r > 2, and M G L"^(f2)^^^, with 
m > 1. Let <^ be a divergence free function in C^(f2)^. Putting p as test function in the limit 
Stokes equation (3.8) and using the representation formula (3.4) we have 



/ Du : Dp dx + {Du)^v ■ pdx — / {v <Siu) : Dp dx + Mu ■ p dx 
Jn Jn Jn Jn 

= if, '")//-l(f7)2,Hl(f7)2- 



(3.52) 



Due to the regularity of Q the set of divergence free functions is known to be dense in the space 
of divergence free functions in HQ^flY (see, e.g., [17]). Moreover, by the higher integrability of 
V and M the mapping 



/ {Du)^v ■ pdx — {v i^u) : Dp dx + Mu ■ p dx 
Jn Jn Jn 



is continuous in i^Q (f2)^. Therefore, considering in (3.52) a divergence free strong approximation 
p oi u in HQ{Qy we get 

/ \Duf dx + {Du)^v-udx— I [v ®u) : Dudx + I Mu ■ udx = {f ,u) H-'i^{ny,H'^{nY ■, 
Jn Jn Jn Jn ' ° 

which is (3.21). This equality clearly implies the uniqueness of a solution u G Hq{Q,Y of (3.8), 
with Mu -ue L\n). 

It remains to prove the corrector result (3.22). Let p G C^{Q). Applying successively the 
triangle inequality and the Cauchy-Schwarz inequality we have 



/ 


\Du,- 


Jn 






[ \Du, 




In 




f \Du, 




In 




[ \Du, 




In 



hence by the boundedness of in ifo(^^)^, 

/ \Due- Du- D{Weu)\dx < / \Due- Du- D{Wep)\dx + c\\u- p\\Hi(n)^. (3.53) 
Jn Jn ° 

On the other hand, proceeding as in fourth step owing to the second convergences of (3.20) 
and (3.10) (which hold in the weak-* sense of measures on Q) we get similarly to (3.48) the 
equality 

/ \Due- Du- D{Wep)\'^ dx ^ {f,u)H-i{n)^,H^{n)^ - / \Duf dx 
Jn Jn 



n 



-2 / Mu-pdx+ I M p ■ pdx + o{l). 
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Hence, taking into account equality (3.21) and using the Holder inequality combined with the 
embedding oi Hq(Q) in any L^{Q) space, it follows that 

\Du,- Du- D{W,ip)f dx = M {u - If) ■ {u - (fi) dx + o{l) 
Jn Jn (3.54) 

< c ||M||i™(f^)2x2 \\u - ^\\Hi(n)2 + o(l). 
Therefore, by (3.53) and (3.54) we obtain the inequality 

limsup / \Dus — Du — D{WsU)\dx < c\\u — ip\\ffi(^Q'^2 , for any 99 G Cj?°(r2)^, (3.55) 

which implies the desired convergence (3.22) and concludes the proof of Theorem 3.3. □ 

As in the scalar case we show in the next section that the equi-integrability condition 
is crucial to derive the limit Brinkman equation (3.8) with the matrix-valued function M 
introduced by L. Tartar [14, 16]. 



3.3 A counter-example 

Let Q be a regular bounded domain of R^. For £ > 0, let be the intersection of Q with the 
periodic lattice of disks of center 2e k, k & 7? , and of radius eVe such that 

47r 

7 e (0,00). (3.56) 



£^ I In I £->-o 

This geometry was used by Cioranescu, Murat [5] for the Laplace equation and by Allaire [2] for 
the Stokes equation, in order to derive a "strange term" of zero-order from the homogenization 
of the Dirichlet boundary conditions on the small disks. 

In the square Y := (—1,1)^, let Q be the disk centered at the origin and of radius 1, and 
let Qr^ be the disk of same center and of radius with measure {Qr^l = t^^e- Then, for 
/ G H~^{fl)^, we consider the Stokes equation 

- Aii, + 1^ Jii, + Vpe =/ inQ 

div(M,) =0 inQ (3-57) 

Ue =0 on dQ. 

Note that, in view of the definition of cu^, we have Icu^l ~ \ {Qr^l- Moreover, if G HQ[fl) is 
the solution of the Laplace equation 

A^, = i^ in2)'(Q), (3.58) 

we have 

-^ = curl(vs) in D'(0), where Ve'.^JVze. (3.59) 

Hence, the Stokes problem (3.57) is of the same type as (3.6). On the other hand, using 
successively the Cauchy-Schwarz inequality and the estimate (3.67) below combined with (3.56) 
we have 

/ \V Ze\'^ dx ^ - -f ZEdx<j^ (-f z'^dx] < c || V^e||L2(n)2, 
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which imphes that is bounded in R^i^. Therefore, the sequence is bounded in L^(f2)^. 
Moreover, since by periodicity the sequence converges weakly-* to \ in M(J1), we get 

— ^ f; weakly in (0)2, with curl (^;) = J in (3.60) 

On the other hand, it is not difficult to check that is not equi-integrable in L'^i^'^ . In fact, 
the following result shows that Theorem 3.3 does not hold for this particular sequence v^: 

Theorem 3.5. The sequence weakly converges in ifo(^)^ '^o the solution u of the Brinkman 
equation 

Am + - Jm + Vp + Vu = f in fl 

div{u) =0 inn (3-61) 
u — on dD,, 
where the extra zero-order term Vu is given by 

^'^^ ^ ^ Jue^ Vu weakly-* in M{nY, (3.62) 



and r is the constant matrix defined by 

Moreover, the matrix obtained from convergence (3.20) according to the Tartar approach is 
given by 

M = ^I. (3.64) 
47 

Remark 3.6. The matrix F of the Brinkman equation (3.61) is not symmetric contrary to the 
matrix M arising in the Tartar approach. Moreover, we have 

Fu-u < Mu -u if M 7^ 0. 

The gap between the two previous energies (which are the energies dissipated by viscosity 
according to [13]) is due to the loss of equi-integrability of the sequence defined by (3.59). 
Therefore, the equi-integrability of can be regarded as an optimal condition to ensure the 
result of Theorem 3.3. 

Remark 3.7. It is worth to mention that the pathology displayed in Theorem 3.5 is not due 
to the absence of correctors. Indeed, with the oscillating sequences vl, vl defined by (3.68), 
(3.69) below, the following corrector result holds: 

Proposition 3.8. Assume that u e W^''^{Q)'^ for some r > 2. Then, we have 



Ue — u — vivl — V2 vl — > strongly in H^{fl), 

\_, . . (3-65) 

7^ + 1 



where v = (vi, V2) := ^ , ^ (- Ui + 71^2, -U2- ^Ui) . 



Remark 3.9. If the right-hand side / belongs to W^'^''^{Qy for some r > 2, then using the reg- 
ularity results for the Stokes equation (see, e.g., [9]) the solution u of the Stokes equation (3.61) 
belongs to W^''^{0.y. This provides a quite general condition under which the strong conver- 
gence (3.65) holds. 
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The proof of Theorem 3.5 is partially based on the properties of the test functions vl, 
defined by (3.68), (3.69) below, and introduced by Allaire [2]. They were also used in [4] to 
derive a homogenized Brinkman type equation but, contrary to (3.6), from a Stokes equation 
without zero-order term. More precisely, in [2] the velocity is assumed to be zero in the set u^. 
In [4] the viscosity is assumed to be very high in cylinders of section cUe, which leads to a three- 
dimensional nonlocal Brinkman equation. In the perturbed Stokes equation (3.57) a highly 
oscillating zero-order term is concentrated on w^. 

On the one hand, the sets Qr^ and cu^ satisfy the following estimates: 
Lemma 3.10. There exists a constant C > such that 



\/V e H\Y), 



/ Vdy--f 



Vdy 



<C^y\lnr,\\\VV\\L2 



J We 



(3.66) 



(3.67) 



Proof. Estimate (3.66) can be easily proved using the polar coordinates. Estimate (3.67) is an 
immediate consequence of the Lemma 3 of [11], and can also be deduced from (3.66). □ 

On the other hand, consider the eF-periodic functions vl and , for i = 1, 2, defined by 



X 



1 



vl(x):^V:i-), pI(x):=-P:1-), for x G 



X 



where e H\{Y) are PI e L'^{Y) are the F-periodic functions defined by 



Si in Qr 



in/\Q, n^ = o mg..u(r\g), I p:dy = o, 



which solve the Stokes equation 

-AV: + VP:^0 mQr,\Q. 
Moreover, the sequences and P^ satisfy the following estimates: 
Lemma 3.11. There exists a constant C > such that 



(3.68) 



(3.69) 



(3.70) 



|V7lU2(y)2 + ||-DKllL2(y)2x2 



+ \\P. 



£|IL2(Y) 



< 



c 



In I 



(3.71) 



and for any function V e H^{Y), 



< 



DV: :DVdy- J P^ div (V) dy - Ye e» 

C 



V 



Qr 



V 



Y\Q 



(3.72) 



In 



\DV\ 



L2(y)2x2, 



where -f denotes the average value and 





i^*-" 

£->0 In Te 



(3.73) 
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Proof. Estimate (3.71) can be proved using the polar coordinates (see also [2]). Estimate (3.72) 
is a straightforward consequence of the Lemma 3.3 of [4] (with a refinement for the right-hand 
side of the inequality) . □ 

Proof of Theorem 3.5. The proof is divided into two steps. In the first step we determine 
the homogenized Brinkman equation (3.61). The second step is devoted to the computation of 
the matrix M defined in the Tartar approach. 

First step: Determination of the homogenized equation. 

Using as test function we have 

/ l-Dwep = (/,M£)i/-i(n)2,m(f2)2 ^ c ||/||ij-i(f2)2 ||-DM£||L2(n)2x2 , 
Jn 

which imphes that Ue is bounded in Hq{0.)'^. On the other hand, let (p e C^{0.) with zero 
Q-average. There exists (see, e.g., [3]) a vector-valued function $ e C^(Q)^ such that 

div ($) = (p inn and ||$||/ji(n)2 < c ||92||L2(n) , 

where the constant c is independent of Lp, $. Using $ as test function in equation (3.57) and 
applying successively the Cauchy-Schwarz inequality, estimates (3.67) and (3.56) we get 



/ Ps^dx 


< 


Jn 






< ( 




< ( 



if, ^)H-i{n)2,H^(n)2 



n \Qre 



< C ||£)$||i2(Q)2x2 + C 



< C ||i:)$||i2(n)2x2 C£^ I Inr^l ||£)W£||i,2(f2)2x2 ||i:)$||i2(f2)2x2 < C ||(/7||L2(n). 

This combined with the regularity of Q implies that p^ is bounded in L^(f2)/M. Therefore, up 
to a subsequence the following convergences hold 



Ue — ^ u weakly in Hq{Q)^ 
Ps — ^ p weakly in L^(0)/]R. 



(3.74) 



Now, we have to determine the limit of the sequence Ju^. On the one hand, re-scaling 
inequality (3.72) we obtain that the functions vl and pi, i — 1,2, of (3.68) and any function 
V e Hq{0,)'^ satisfy the inequality 



< 



Dvl -.Dv- I pi div (v) - ^ e. 



\Dv\ 



Jn\Qre\ Jn\y\Q\^eJ 



(3.75) 



L2(0)2X2. 



£ I In I 

Moreover, by (3.71) and (3.56) the following convergences hold 

vl — ^ weakly in H\n)^ 
Pi 



weakly in L2(0)/M. 



(3.76) 



Then, applying inequality (3.75) with v — (pu^, (p & C^(Jl), we deduce from (3.73) and (3.56) 
that 

I. 



f Dvi:Du,ip-{^ + o{l))ei 
Jn 



■-We 



n \Qre 



ipu, 



In^ 



ipu 



0(1). 



(3.77) 
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On the other hand, putting ipvl as test function in (3.57), using that vl = Cj in Us and the 
convergences (3.76), (3.74), we have 

/ Du, ■.Dviip+ I Ju, • (/? = o(l). (3.78) 



Denote 



v := hm 777^ weakly-* in M(f2)^, 

£-)-0 



where the hmit holds up to a subsequence by virtue of the estimate (3.67) combined with the 
Cauchy-Schwarz inequality. Then, equating (3.77) and (3.78) and passing to the limit we get 

for i = 1,2, 

/ (pei ■ u — ^ / (p ei ■ Jv + / ipei ■ u, for any (p e C^{Q), 
Jn Jn 4 Jq 

which implies the equality v — ^ Jv + -^u. Hence, we deduce the convergence 

1^ Jii, ^ = I (/ - 7 J)-i u = 4^:^^^ {I + lJ)u weakly-* in M(Q)l (3.79) 

Therefore, passing to the limit in (3.57) with (3.79) we obtain the homogenized equation 

_Ax^+ljx^ + Vp+ ^^^/^^^ {^I-J)u = f iTiD'iSl), (3.80) 

which yields the desired Brinkman equation (3.61) with the matrix F of (3.63). 
Second step: Derivation of the matrix M. 

Let A e R^. Consider the solutions W^^^ e ^10^) (^^® F-periodic functions in 

Hl^^{m?)) and QL G -L?(y)/R of the perturbed Stokes problem 



-AW^^. + ^(^-i)^A + VQt =0 inR^ 

div(iy„y =0 inR2 



Wu\ is y-periodic 



Y 



Note that the first equation of (3.81) is equivalent to the variational formulation in the torus, 

'iVeHliY), J DW^^^,: DVdy + s (J V - j ■ J\- J Ql,div{V)dy^O. (3.82) 

Hence, the re-scaled functions w^^^ and q^^^ defined by 

wl{x):^sWl,l^^) and q^x) -.^ (^^) , for x G O, (3.83) 
are sF-periodic solutions of the problem 

-A<+(^(f)-i)7A + V,t =0 in.. pg^j 

div (ty^^) =0 in 
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First of all, let us determine a priori estimates satisfied by the sequences VFA, wi'^, QL, 



and q^^^. Putting W^^^ as test function in equation (3.82) we have 

/ \DWlXdy + e-f JX-Wl,dy^O, 
Jy Jq,^ 

hence by the estimates (3.66) of Lemma 3.10 and (3.56) 



(3.85) 



|-^^I1,£|L2(Y)2X2 ~ ^ 



Therefore, W^^ is bounded in Hj{Yy, and there exists a constant vector e such that 
up to a subsequence we have 

(3.86) 



On the other hand, let (p G C^{Y) with zero F-average. There exists $ G Cj°°(y)^ with zero 
y-average such that 

div {^) — if in and ||$||i2(y)2x2 < c ||(^||i:,2(y), 

where c is a constant independent of (/?, Putting $ as test function in (3.82) we have by 
(3.66) and (3.56) 



< 



DWt-.D^dy 



Y 



$ dy 



< c ||Z?$||i2(y^2x2 + cey I Inr^l ||£'$||x,2(y)2x2 

< c||D$||i2(y)2x2 < C||(^||i2(y), 



hence is bounded in Lj{Y)/'R. Prom the boundedness and the y-periodicity of W^^^ and 
Q^^^ we thus deduce that the sequences w^^^ and g^^ of (3.83) satisfy the convergences 



4e 



weakly in H\Qf 
weakly in L2(0)/M. 



(3.87) 



Now, let us check that the periodic function w^^^ of (3.84) gives the same matrix M (3.20) as 
the function of (3.7) which satisfies a Dirichlet boundary condition. Since M is symmetric, 
this is equivalent to prove that for any A G R^, 



{DwLyVe ■ A - (Dw^Yv, • A ^ mV'iQ), 



(3.88) 



where is defined by (3.59). Let ip G C^{Q). Putting p^w^^^ in the equation (3.7) satisfied by 
and ipw^ in the equation (3.84) satisfied by w^^, and using the convergences (3.87) satisfied 
by w^^^, as well as the similar ones satisfied by w^, q^, we get 



Dw^ : Dwl, p - I curl {v,) Jwl ■ A ^ 



/ Dw^^^ : Dw^ (f- curl (wj Jw^ ■ A (p 
Jn ' Jn 



£->0 



(3.89) 



0. 
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Moreover, by the representation formula (3.11) we have 

curl (ve) JwL - (DwLV Ve — ^ 

^ ' ^ in 'D'iQf. (3.90) 

curl {ve) Jw^ - {Dw^Y Ve — ^ 

Therefore, combining (3.89) and (3.90) we obtain the desired convergence (3.88). 

It remains to determine the matrix M. On the one side putting as test function in 
(3.84) and using the convergences (3.88), (3.20), and on the other side using the £y-periodicity 
of Wj^g (3.83), we get similarly to (3.9) and up to a subsequence 

\DwlJ^ ^ MA • A and \DwIJ\ lim {J- \DW^^J\ dy^ weakly-* in M(l^). (3.91) 

This combined with (3.85) and (3.86) gives 

MA • A = ^ JW^ ■ A. (3.92) 

Let us compute the constant vector W^. To this end, putting the divergence free function W^^^ 
in the inequality (3.72) satisfied by V^, i — 1,2, and taking into account the estimates (3.71), 
(3.56) and the boundedness of W^^ in H\Yf, we have 

^ DV: : DWl, dy ^Ye^i- W^, dy^ + o{e) . (3.93) 

Moreover, putting the divergence free function in (3.82) with = in Qr^, we get 

/ DW^\ : DV: dy = -e(-f - I V'^ ■JX^eJei-X + o{e), (3.94) 

since by (3.71) strongly converges to zero in L^(y)^. The estimates (3.93) and (3.94) divided 
by e together with (3.86), (3.73) and (3.56) imply that 

-fei-W^ = Jci-X or equivalently = - -JX. (3.95) 

7 

This combined with (3.92) yields the value (3.20) of the symmetric matrix M. □ 

Proof of Proposition 3.8. Let v = (^1,^2) ^ W^''^'{iiy. Considering the functions u*, i = 1, 2, 
which are defined by (3.68) and satisfy the convergences (3.76), we have 

Ee / I Due — Du — vi Dv\ — V2 Dv^\'^ dx 

= / \Due\'^dx- / \Du\'^dx+ / {vflDvll'^ + v^\Dvl\'^) dx (3.96) 
Jn Jn Jn 

-2 {vi Due : Dvl + V2 Du^ : Dvl) dx + o(l). 

Putting Ue in equation (3.57) and u in equation (3.61) we get 

/ \Due\^dx- I \Du\^dx^ {f,u)H-i{Q.Y,Hl{Q.Y- / \Du\'^ dx + o{l) 
Jn Jn Jn 



Fu ■ u dx + 0(1) — ——J- — TT / \u\"' dx + 0(1] 



n 



(3.97) 



4 (7' + 1) Jn 
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Moreover, putting V^' in estimate (3.72) together with — Cj in Qr^, (3.71), (3.73), (3.56), 
and using the £y-periodicity of Dvl, we get 

iDv'f hm (\ -f iDKTdy ) = ^ weakly-* in M(Q), 

hence since Vi e C(n), 

/ {vj \Dvl\^ + vl \Dvl\') dx = j[ \v\^ dx + o(l). (3.98) 

Estimates (3.96), (3.97) and (3.98) thus imply that 

E^^—^ / \u\'^dx + - / \v?dx-2 / iviDu^: Dvl + V2Du^: Dv't)dx + oil). (3.99) 

4(7'+!)^ 4 in ^ 

On the other hand, applying the estimate (3.75) with the function v — ViU^, i — 1, 2, and using 
the convergences (3.76), (3.73), (3.56) and (3.79), we obtain 



Dvl : Due Vidx = Dvl : D{vi u^) dx + o(l) 
n Jn 



£2 



7^ 



4(72 + 

This combined with (3.99) yields 

T / L.12 , , 7 / , .12 



/ 64 ■ i'y I — J) uvidx — / UiVidx + oil 
l)Jn 4 7n ^ 



)■ 



= wo', 7 / \u\'''dx + - I \v\'''dx 

1) Jn 



^ (3.100) 



4(72 + 1)^' ' 4 

+ ^(^/n"-''" + ^(7TI)//"-^'^ + ^(^)- 
Putting the function 

V := -(/ + 7 J)w 

7"^ + 1 

in (3.100) we get 

E^= \Due- Du-ViDvl-V2Dvl? dx — > 0. (3.101) 
Jo. ^-^-0 

Finally, since the sequences vl strongly converge to zero in (fi)^ by (3.76) and u e W^''^{VtY, 
the strong convergence (3.65) is a straightforward consequence of (3.101). □ 



References 

[1] S. Alinhac & P. Gerard: Pseudo- differential operators and the Nash-Moser theorem, 
translated from the 1991 French original by Stephen S. Wilson, Graduate Studies in Math- 
ematics 82, American Mathematical Society, Providence, RI, 2007. 



30 



[2] G. Allaire: "Homogenization of the Navier-Stokes equations in open sets perforated 
with tiny holes I. Abstract framework, a Volume distribution of holes", Arch. Rational 
Mech. Anal, 113 (1991), 209-259. 

[3] M.E. BOGOVSKL "Solution of the of the first boundary value problem for the equation 
of continuity of an incompressible medium", Soviet Math. Dokl., 20 (1979), 1094-1098. 

[4] H.C. Brinkman: "A calculation of the viscous force exerted by a flowing fluid on a dense 
swarm of particles, Appl. Sci. Res., Al (1947), 27-34. 

[5] D. ClORANESCU & F. Murat: "Un terme etrange venu d'ailleurs, I & 11", in Nonlinear 
partial differential equations and their applications, College de France Seminar, II &c III, 
ed. by H. Brezis & J.-L. Lions, Research Notes in Math., 60 &c 70, Pitman, London, 
1982, 98-138 & 154-178. English translation: "A strange term coming from nowhere", 
in Topics in the mathematical modelling of composite materials, ed. by A. Cherkaev & 
R.V. Kohn., Progress in Nonlinear Differential Equations and their Applications, 31, 
Birkhauser, Boston 1997, 44-93. 

[6] G. Dal Maso & A. Garroni: "New results on the asymptotic behaviour of Dirichlet 
problems in perforated domains". Math. Mod. Meth. Appl. Sci., 3 (1994), 373-407. 

[7] P. Gerard: "Micro-local defect measures" Com. Part. Diff. Equ., 16 (11) (1991), 1761- 
1794. 

[8] G. Geymonat : "G. Sul problema di Dirichlet per le equazioni lineari ellittiche." , Ann. 
Scuola Norm. Sup. Pisa, 16 (3) (1962), 225-284. 

[9] O.A. Ladyzhenskaya: The mathematical theory of viscous incompressible flow, trans- 
lated from the Russian, Mathematics and its Applications, 2, Gordon and Breach, Science 
Publishers, New York- London- Paris 1969. 

[10] F. Murat: "Compacite par compensation", Ann. Scuola. Norm. Sup. Pisa, Serie IV, 5 
(3) (1978), 489-507. 

[11] C. PiDERl & P. Seppecher: "A second gradient material resulting from the homoge- 
nization of an heterogeneous linear elastic medium". Continuum Mech. and Thermodyn., 
9 (5) (1997), 241-257. 

[12] E. Sanchez-Palencia: Non Homogeneous Materials and Vibration Theory, Monographs 
in Physics, 127, Springer- Verlag Berhn 1980. 

[13] L. Tartar: "Homogeneisation en hydrodynamique" , in Singular Perturbation and Bound- 
ary Layer Theory, Lecture Notes in Mathematics, 597, Springer, Berlin-Heidelberg 1977, 
474-481. 

[14] L. Tartar: "Remarks on homogenization", in Homogenization and effective moduli of 
materials and media, IMA Vol. Math. Appl., 1, Springer, New- York 1986, 228-246. 

[15] L. Tartar: "H-measures, a new approach for studying homogenization, oscillations and 
concentration effects in partial differential equations", Proc. Roy. Soc. Edinburgh 115 A 
(1990), 193-230. 

[16] L. Tartar: The General Theory of Homogenization: A Personalized Introduction, Lecture 
Notes of the Unione Matematica Italiana, Springer- Verlag, Berhn Heidelberg 2009, pp. 471. 



31 



L. Tartar: "Topics in nonlinear analysis", Publications Mathematiques d'Orsay, 78 (13), 
Orsay 1978, pp. 271. 



32 



